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We proposed two different methods for simulating theta-vacuum systems. Both can be affected by systematic 
effects but these errors can be controlled by comparing the results of the two unrelated methods. This has been 
done for some analytically soluble models and here the same procedure is successfully applied to the interesting 
case of CP 9 in the scaling region. 



Quantum Field Theory with topological terms 
is a special case of complex Euclidean actions: 
the action is naturally divided in a real part and 
a pure imaginary one, the latter in the form 
of i times 9 times an integer that depends on 
the gauge configuration. The same structure 
is shared by spin systems coupled to imaginary 
magnetic field: they do not describe physical sys- 
tems but are nevertheless useful to improve our 
understanding of complex actions. 

For this class of models we proposed originally 
a method for reconstructing the 9 dependence 
from real action simulations. The main ingredi- 
ent was to compute the probability distribution 
function of the topological charge at 9 = from 
simulations at imaginary 9 — ih 1 . An ana- 
lytical form of the p.d.f. is extracted from the 
measurement of the mean value of the density of 
topological charge as a function of h. This is done 
using the saddle point approximation and doing a 
fit of numerical data with a suitably chosen func- 
tion. At the end a multi-precision code is used 
to compute the sum over the different topological 
sectors and calculate the relevant physical quan- 
tities. The arbitrariness of the fitting function as 
well as the use of saddle point relation (valid up to 
0{\/V) corrections at finite volume) is a possible 
source of uncontrolled systematic effects. This is 
not the case in some soluble models considered 
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(2D compact U(l) at any coupling as well as ID 
AF Ising models coupled to imaginary magnetic 
field) where the results where reproduced at % 
level (or better) but in general it is not possi- 
ble to put any upper bound on the magnitude of 
these systematic effects. 

This was precisely our motivation to construct 
an alternative and independent approach [2]. 
This second approach requires as input the same 
data as the first one i.e. x{h) 1 the measurement 
of the mean value of the density of topological 
charge as a function of h. This time we do not 
try to reconstruct the p.d.f. but use the fact that 
the sum over topological sectors that defines the 
partition function can naturally be written as an 
even polynomial in the variable z = cos | and 
the derivative of the free energy density respect 
to logz 2 gives y(z) = x(0)/tan| where x(9) is 
the topological charge density as function of 9. 
From the last relation and its analytical contin- 
uation to imaginary 9 we see that the behavior 
of y(z) for z — cosh(/i/2) S [1, oo) is related to 
x(h) at h > while for z = cos(0/2) 6 [0, 1] 
we can get x(9). From real action simulations 
we can get y(z) for any value of z > 1, i.e. 
up to y(z = 1) = 2xo = ^ dx ^P I e=o an( i since 
y(z -> 0) = lime^ ta ^ 2) the behavior of y(z) 
near the origin is related to x(9 ~ it). Using 
the transformation y\{z) = y(e x l 2 z) we can de- 
fine an effective exponent 7A = f log t n£ d 
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Figure 1. The effective exponent 7a. Lower curve: 
the symbols are from simulation and the continu- 
ous line is the quadratic fit of data. Upper curve: 
the continuous line is the analytic result for imag- 
inary 9 and the dotted part corresponds to real 9 
solution. 



can help to distinguish between different behav- 
iors at 9 ~ 7r since if q(9 — > 7r) oc (it — 6) a then 
7a (y 0) = 1 + a. If 7 A(y -> 0) = 1 CP sym- 
metry is spontaneously broken at 9 — n, values 
between 1 and 2 indicate a CP symmetric vac- 
uum with a divergent susceptibility (second order 
phase transition) and j\ (0) = 2 analiticity of the 
free energy at 9 = n. It is possible to measure 
7a only for values of y larger than 2\o and our 
experience with several models is that this is a 
smooth function of the argument. 

Our second method is based on the feasibil- 
ity of a simple extrapolation procedure inside the 
forbidden region to reconstruct the full 7a depen- 
dence. Once we have the effective exponent at 
any y > we can reconstruct the topological 
charge density dependence on 9 using an itera- 
tive procedure as explained in [2|. 

As an example we plot in Fig. 1 the effective 
exponent for ID AF Ising model (CP symme- 
try is spontaneously broken at 9 = it) and for 
a toy model that resembles the free instanton 



gas solution (the partition function is Zy(9) = 
(1 + ^Icos^)^ where A is a constant and the CP 
symmetry is enforced at 9 = 7r). In both cases 
a simple quadratic fit allows to extract values of 
7a (y = 0) compatible with 1 and 2 respectively 
(with error at .1% level). Also the reconstruction 
of the order parameter gives results with errors 
at level of few % maximum (see [2])- 

This shows that the extrapolation of the ef- 
fective exponent is a reasonable possibility. Of 
course in doing that we assume a smooth behav- 
ior near the origin and this is equivalent to as- 
sume that no phase transitions are present as a 
function of 9 except at most at 9 = it. We can 
also naively expect that the extrapolation is safer 
for asymptotically free models since in this case 
the numeric data can extend to smaller values of y 
since the topological charge susceptibility \o goes 
to zero increasing 0. Since A is a free parameter 
in the game an internal consistency check of the 
procedure is to verify the independence of physi- 
cal observables on it. 

Clearly, once again, systematic effects intro- 
duced by the arbitrary choice of the extrapolat- 
ing function can be present. What is important is 
that we expect the systematics to be different for 
the two different methods. On the other hand, if 
both methods give consistent results, our confi- 
dence on the quality of the results is very high. 

Once we have two methods to crosscheck our 
results the main goal would clearly be to under- 
stand the role of the 9 parameter in QCD and 
its connection with the strong CP problem. The 
application of our methods to QCD with theta- 
tcrm would be straightforward provided we have 
a practical definition of the topological charge in 
order to do simulations at imaginary 9. It turns 
out not to be the case: the known regularizations 
of the topological charge (e.g. ginsparg-wilson) 
are so involved that at the moment they are not 
of practical use. A good possibility for studying 
the #-term effect on the low energy dynamics are 
two dimensional CP^ -1 models since they have 9- 
vacua, a well defined topological charge and many 
features in common with four-dimensional Yang- 
Mills theories: they are asymptotically free, pos- 
sess instanton solutions, and admit a 1/N expan- 
sion. 
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Figure 2. The topological charge density in L = 
200 lattice. The points refer to the first method, 
the diamonds to the second one (statistical errors 
are smaller than the symbol). 



We already considered CP 3 in the strong cou- 
pling region and are now interested to extend 
these studies to the continuum limit. In this case 
the effect of dislocations has to be negligible not 
to spoil the scaling properties of the theory. From 
6 = simulations we know that N — 10 is large 
enough [3] thus we concentrated our efforts on the 
scaling region of the CP 9 model. 

We have adopted for the action the "auxiliary 
U(l) field" formulation and the topological charge 
operator is defined directly from the U(l) field 
(see 0). For N = 10 and h = the scaling win- 
dow starts at j3 around 0.75 and we decided to 
concentrate our efforts at (3 =0.8 , 0.85 , 0.90. 
The lattice sizes considered are L = 100, 200 for 
the two larger couplings and L = 200 for the 
smaller one. A standard Metropolis algorithm 
was used to generate O(10 6 ) configurations for 
each value of (3 and h (we refer to for a full 
discussion of those details). 

In Fig. 2 we compare the results of the two 
methods for two values of (3: the agreement is at 
(few) % level, compatible with the statistical one 
(the errors in the plot are the statistical ones). 



Figure 3. The topological charge density divided 
by the measured topological susceptibility at 8 = 
for the L = 200 lattice. 



Following our previous consideration this can be 
considered the strongest possible a posteriori con- 
sistency check of the results. 

We used our high quality data to address the 
scaling properties of the theory too. This has 
been done perturbatively and non perturbatively 
@|. In Fig. 3 we show how the non perturba- 
tive scaling (only data from the first method are 
plotted) is realized up to violations that are at 
most ±4%. As expected the biggest violations 
are for the smaller inverse coupling. We can con- 
clude that we have strong indications that, for 
this model, CP is spontaneously broken in the 
continuum, as predicted by the large N expan- 
sion. 
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